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We investigate the conductance spectra of a normal/superconductor graphene junction using the extended 
Blonder-Tinkham-Klapwijk formalism, considering pairing potentials that are both conventional (isotropic s- 
wave) and unconventional (anisotropic d-wave). In particular, we study the full crossover from normal to specu- 
lar Andreev reflection without restricting ourselves to special limits and approximations, thus expanding results 
obtained in previous work. In addition, we investigate in detail how the conductance spectra are affected if it is 
possible to induce an unconventional pairing symmetry in graphene, for instance a d-wave order parameter. We 
also discuss the recently reported conductance-oscillations that take place in normal/superconductor graphene 
junctions, providing both analytical and numerical results. 
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I. INTRODUCTION 

A key issue in understanding low-energy quantum transport 
at the interface of a non-superconducting and superconduct- 
ing material, e.g. a normal/superconductor (N/S) interface, is 
the process of Andreev reflection. Although the existence of 
a gap in the energy spectrum of a superconductor implies that 
no quasiparticle states may persist inside the superconductor 
for energies below that gap, physical transport of charge and 
spin is still possible at a N/S interface in this energy-regime 
if the incoming electron is reflected as a hole with opposite 
charge. The remaining charge is transferred to the supercon- 
ductor in the form of a Cooper pair at Fermi level. The study 
of Andreev reflection and its signatures in experimentally 
observable quantities such as single-particle tunneling and 
the Josephson current has a long history (see e.g. Ref. [l|and 
references therein). Only recently, however, has this field of 
research been the subject of investigation in graphene N/S 
interfaces^i^. 

Graphene is a monoatomic layer of graphite with a hon- 
eycomb lattice structure, as shown in Fig. [1] and its recent 
experimental fabrication^-'' has triggered a huge response in 
both the theoretical and experimental community over the last 
two years. The electronic properties of graphene display sev- 
eral intriguing features, such as the six-point Fermi surface 
and a Dirac-like energy dispersion, effectively leading to an 
energy-independent velocity and zero effective mass at the 
Fermi level. This obviously attracts the interest of the theo- 
rist, but graphene may also hold potential for technological 
applications due to its unique combination of a very robust 
carbon-based structural texture and its peculiar electronic fea- 
tures. 

Condensed matter systems with such 'relativistic' elec- 
tronic structure properties as graphene constitute fascinat- 
ing examples of low-energy emergent symmetries; in this 
case, Lorentz-invariance. At half-filling, the Fermi level of 
graphene is exactly zero which renders the Fermi surface to 
be reduced to a six single points due to the linear intersection 
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FIG. 1: (Color online) a) Sketch of real-space lattice of graphene, 
consisting of two hexagonal sublattices A and B. The interatomic dis- 
tance a is equal between all lattice points, b) k-space (momentum- 
space) lattice of graphene, including the hexagonal Brillouin zone. 
Only two inequivalent points exist on the BZ boundary, termed K 
and A". 

of the energy bands (see Fig. [3]and|4|i. The linear disper- 
sion relation is a decent approximation even for Fermi levels 
as high as 1 eV, such that the fermions in graphene behave 
like they are massless in the low-energy regime. The fact that 
the fermions around Fermi level obey a Dirac-like equation 
at half-filling introduces Lorentz-invariance as an emergent 
symmetry in the low-energy sector. Another example where 
Lorentz-invariance appears for low-energy excitations is in 
one-dimensional interacting fermion systems, where phenom- 
ena like breakdown of Fermi-liquid theory and spin-charge 
separation take place. When Lorentz-invariance emerges in 
the low-energy sector of higher-dimensional condensed mat- 
ter systems, it is bound to attract much interest from a funda- 
mental physics point of view. Another interesting feature of 
graphene are the nodal fermions that are present at the Fermi 
level at half-filling. When moving away from half-filling by 
doping, the excitations at the Fermi level are no longer nodal. 
The nodal fermions of graphene hold certain similarities to, 
but also important differences from, the nodal Dirac fermions 
appearing in the low-energy sector of the pseudogap phase of 
d-wave superconductors such as the high-Tc cuprates. In con- 
trast to graphene, the nodal fermions in the high-Tc cuprates 
track the Fermi level when these systems are doped and thus 
represent a more robust feature than in graphene. For an illus- 
tration of the latter scenario, consider Fig. |2] which contains 
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a sketch of the Fermi surface in the cuprates when including 
terms up to next-nearest neighbor hopping. The nodal lines 
of the (i^2_y2-gap intersect the Fermi surface at exactly four 
points, which permits the existence of nodal fermions at those 
points in k-space. However, in contrast to graphene, doping 
the system will in this case simply move the position of the 
Fermi arc with respect to the nodal line of the superconduct- 
ing gap, such that the nodal fermions persist in the system^. 
Nonetheless, the existence of the Dirac cones in graphene rep- 
resents an important example of emergent non-trivial symme- 
tries at long distances and low energies in higher (more than 
one) dimensional systems. 



Fermi surface Nodal fermions 




FIG. 2: (Color online) Sketch of the Fermi surface and anisotropic 
d-wave gap as believed to be present in high-Tc cuprate supercon- 
ductors. The nodal fermions reside at the intersection of the nodal 
lines of the gap and the Fermi surface. 

Although superconductivity does not appear intrinsically 
in graphene, it may be induced by means of the prox- 
imity effect by placing a superconducting metal electrode 
near a graphene layer^''''^'"''^. Recent theoretical work^ '^ 
have considered coherent quantum transport in N/S and nor- 
mal/insulator/superconductor (N/I/S) graphene junctions in 
the case where the pamng potential is isotropic, i.e. s-wave 
superconductivity. However, the hexagonal symmetry of the 
graphene lattice permits, in principle, for unconventional or- 
der parameters such as p-wave or d-wave (characterized by a 
non-zero angular momentum of the Cooper pair). A complete 
classification of the possible pairing symmetries on a hexago- 
nal lattice up to /-wave pairing {I = 3) was given by Mazin 
and Johannes'"*, with the result given in Tab. |T] We underline 
that the notation "insulator" in this context refers to a normal 
segment of graphene in which one experimentally induces an 
effective potential barrier Vq. As we shall see, such a potential 
barrier has dramatically different impact upon the transport 
properties in graphene as compared to the metallic counter- 
part. 

The intrinsic spin-orbit coupling in graphene is very weak, 
as dictated by the low value of the carbon atomic number, such 
that we will neglect it in this work. We will also disregard the 
electrostatic repulsion as mediated by the vector potential A. 
At first sight, this might seem as an unphysical oversimpli- 
fication since there is no metallic screening of the Coulomb 



TABLE I: List of all superconducting pairing states allowed for a 
hexagonal lattice up to d-wave pairing, adapted from Ref. I^. An 
orbital angular momentum quantum number / — 0, 1, 2 is denoted 
s,p, d-wave, respectively. For the triplet states (p-wave), the order 
parameter has multiple components, and is conveniently represented 
as a vector dk. 
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interaction in graphene. In an ordinary metal, the renormal- 
ized Coulomb-potential reads V{r) — Vb(?')c^'"/'^, where 
A EE [N{Ef)]-^l'^ is the Thomas-Fer mi screening length and 
N{Ep) is the density of states (DOS) at Fermi level. Since 
pure graphene has zero DOS at Fermi level, one might quite 
reasonably suspect that the screening of charge vanishes, and 
it might seem paradoxical that Coulomb interactions can be 
neglected. The resolution to this is found by realizing that one 
may disregard the Coulomb interaction if it is weak compared 
to the kinetic energy in the problem. Due to the linear disper- 
sion, the kinetic energy is governed by the Fermi velocity vp 
which formally diverges near Fermi level. The divergence is 
logarithmic and precisely due to the Coulomb interactioni^J^. 
The limiting velocity in graphene (due to e.g. Umklapp pro- 
cesses) is nevertheless of order 0(10^) m/s, see e.g. Ref. [TtI 
This is roughly 100 times larger than in a normal metal, and 
it is thus safe to neglect the Coulomb interaction compared 
to the kinetic energy in graphene. In graphene, the Coulomb 
interaction self-destructs. 

In this work, we will study in detail how an anisotropic or- 
der parameter induced in graphene will affect quantum trans- 
port in a N/S and N/I/S junction, extending the result of 
Ref. [iSi In equivalent metallic junctions, it is well-known— 
that the zero bias conductance peak (ZBCP) is an experimen- 
tal signature of anisotropic superconductivity in clean super- 
conductors with nodes in the gap. This is a consequence of 
bound surface states with zero energy at the interface that form 
due to a constructive phase-interference between electron- 
like and hole-like transmissions into the superconductor^^. In 
graphene junctions with superconductors, as we shall see, a 
new phenomenology comes into play with regard to the scat- 
tering processes that take place at the N/S interface. It is 
therefore desirable to clarify how anisotropic superconductiv- 
ity is manifested in the conductance spectra of such a junc- 
tion, and in particular if the same condition for formation of a 
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ZBCP holds for graphene junctions as well. As first shown in 
Ref.JS^ we will demonstrate that in N/I/S graphene junctions, 
novel conductance-oscillations as a function of bias voltage 
are present both for s-wave and d-wave symmetry of the su- 
perconducting condensate due to the presence of low-energy 
'relativistic' nodal fermions on the N-side. The period of the 
oscillations decreases with increasing width w of the insu- 
lating region, and persists even if the Fermi energy in I is 
strongly shifted. This contrasts sharply with metallic N/I/S 
junctions, where the presence of a potential barrier causes 
the transmittance of the junction to go to zero with increas- 
ing w. The feature of conductance-oscillations is thus unique 
to N/I/S junctions with low-energy Dirac-fermion excitations. 
Moreover, we contrast the N/S or N/I/S conductance spectra 
for the cases where s-wave and (ia,2_j,2-wave superconductor 
constitutes the S-side. The former has no nodes in the gap and 
lacks Andreev bound states. The latter has line-nodes that al- 
ways cross the Fermi surface in the gap, and thus features in 
addition to Andreev bound states, also nodal relativistic low- 
energy Dirac fermions. The quantum transport properties in 
a heterostructure of two such widely disparate systems, both 
featuring a particular intriguing emergent low-energy symme- 
try, is of considerable importance. 

This paper is organized as follows. In Sec. [Ill we establish 
the theoretical framework which we shall adopt in our treat- 
ment of the N/S graphene junction. The results are given in 
Sec. [in] where we in particular treat the role of the barrier 
strength and doping with respect to how the conductance is 
influenced by these quantities. In addition, we investigate the 
role of a possible unconventional pairing symmetry induced 
in graphene. A discussion of our findings is given in Sec. FVl 
and we summarize in Sec. |Vl] We will use for 4x4 ma- 
trices, and .?. for 2x2 matrices, with boldface notation for 
three-dimensional row vectors. 



II. THEORETICAL FORMULATION 

A. General considerations 

The Brillouin zone of graphene is hexagonal and the en- 
ergy bands touch the Fermi level at the edges of this zone, 
amounting to six discrete points. Out of these, only two are 
inequivalent, which are conventionally dubbed K and K' , and 
referred to as Dirac points. The band dispersion of graphene 
was first calculated by Wallace^l, and reads 

E = ±7o[l+4cos (^) cos (^)+4cos2 (^)] 

(1) 

where 70 ~ 2.5 eV, and the ± sign refers to the anti- 
bonding/bonding TT-orbital. The remaining three valence elec- 
trons are in hybridized sp^ a-bonds. The energy dispersion in 
the Brillouin zone is plotted in Fig. |3] which reveals the con- 
ical structure of the conduction and valence bands at the six 
Fermi points. The cosine-like conduction and valence bands 
are made up by a mixture of the energy bands from the A- and 
B-sublattices in graphene (Fig. [U, which are linear near the 



Fermi level. This gives rise to the conical energy dispersion at 
the Dirac points K and K' . 




FIG. 3: (Color online) a) The energy dispersion for graphene in the 
Brillouin zone. The upper band is the antibonding 7r-orbital, while 
the lower band is the bonding 7r-orbital. It is seen that the bands touch 
at Fermi level (_Ef = 0) at six discrete points, which constitutes the 
effective Fermi surface, b) Contour plot of the dispersion relation, 
clearly showing the hexagonal structure of the Fermi points. The 
center of each red drop-like structure represents either K 01 K' . 

In order to introduce the new phenomenology of scattering 
processes in N/S graphene junctions, it is instructive to com- 
pare it with the metallic N/S junction. This is done in Fig. |4] 
In the metallic case, an incident electron with energy E < A 
measured from the Fermi energy Ep can not transmitted into 
the superconductor since there are no available quasiparticle 
states. Instead, it is reflected as a hole, represented as a quasi- 
particle with energy E in the hole-like band, such that the 
leftover charge 2e is transferred into the superconductor as a 
Cooper pair at Fermi level. The hole has negative mass, en- 
ergy, wave-vector, and charge compared to the electron which 
is absent. Strictly speaking, only at E = are the wave- 
vectors exactly related through fee — kh since one in general 
has 

fee = y/2m{EF + E), K = y/2m{EF-E). (2) 

At finite energies, the electron-hole coherence will therefore 
be lost after the hole has propagated a distance L ~ 1/E. At 
energies E > A above the gap, Andreev reflection is severly 
suppressed since direct tunneling into quasiparticle states is 
now possible. 

In graphene, a new phenomenology of Andreev reflection 
is at hand due to the band structure which effectively looks 
like that of a zero-gap semiconductor (see also Ref. 3). Since 
the conduction and valence bands touch at the Fermi energy 
Ep = 0, one may distinguish between three important cases: 
i) undoped graphene with Ep — 0, ii) doped graphene with 
Ep > 0, and Hi) heavily doped graphene with Ep ^ 0. These 
different scenarios are shown in Fig. |4] 

In undoped graphene, with Ep ^ A, an incident electron 
with energy E is denied access as a quasiparticle into the 
superconductor, and physical transport across the junction is 
thus manifested through reflection as a hole. When Andreev 
reflection takes place, the transmitted Cooper pair is located 
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at the Fermi level of the superconductor. Energy conservation 
then demands that the missing electron in the normal region 
that is reflected as a hole must be located at —E due to the 
energy conservation, i.e. in the valence band. This is different 
from normal Andreev reflection, since in that case both the 
electron and hole belong to the same band (conduction). For 
specular Andreev reflection, however, they belong to different 
bands. The use of the term "specular" in order to characterize 
this type of Andreev reflection originates with the fact that the 
group velocity Vg and momentum k have the same sign for a 
valence band hole, while in contrast Vg and k have opposite 
signs for a conduction band hole. To see this, consider first a 
usual metallic parabolic dispersion E = k^/2m — Ep for the 
electrons, such that one readily infers from Vg = Vki? that 
Vg = k/m. Therefore, a hole created at a given energy E will 
have Vg = — k/m, since holes have opposite group velocities 
of the electrons for a given wave-vector k. For normal 
Andreev reflection, the holes are located in the conduction 
band and therefore satisfy Vg 1 1 — k. 



In the case of specular Andreev reflection for undoped 
graphene {Ep = 0), a hole is generated in the valence band. 
Since in the valence band, the electronic dispersion reads 
E = — uplkj, the group velocity of electrons is opposite 
to their momentum. Conversely, the group velocity for 
valence holes is parallell to their momentum. This is the 
mechanism behind specular Andreev reflection. In doped 
graphene (Ep > 0), the Andreev reflection can be normal or 
specular, depending on the energy of the incoming electron, 
as sketched in Fig. |4] In heavily doped graphene, [Ep ^ A), 
only normal Andreev reflection is present for subgap energies 
since the distance from Fermi level to the valence band is too 
large for specular AR to occur In the regime Ep G [0,A], one 
has either normal or specular Andreev reflection, depending 
on the incident electron energy E. 



We also comment on the effect of Fermi vector mismatch 
(FVM). Blonder and Tinkham.^^ showed that in a metallic 
N/S junction, a FVM would act as a source for normal re- 
flection, such that one could effectively account for it simply 
by choosing a higher value for the barrier strength Z. Inter- 
estingly, in a ferromagnet/superconductor junction the effect 
of a FVM could not be reproduced by simply shifting Z to a 
higher value, as discussed by Zutic and Valls^^. In the absence 
of an exchange energy, however, the effect of FVM can be 
thought of as a reduction of the Fermi surface that participates 
in the scattering processes, as illustrated in Fig. |5] One may 
parametrize the FVM by the parameter k = kp /qp where kp 
(qp) is the Fermi momentum in the normal (superconducting) 
part of the system. In particular, it is seen that for k > 1, there 
is only possible transmission of quasiparticles (although these 
decay exponentially) up to a critical angle less than 7r/2^. 

Having established the states that participate in the scatter- 
ing at the interface, we now turn to equations that describe 
these quasiparticle states. 
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FIG. 4: (Color online) Graphical illustration of the different scatter- 
ing processes that may take place at a /) metallic and ii) graphene 
N/S junction. While the Andreev reflected hole in case ;) retraces the 
trajectory of the incoming electron, the hole in case ;/) may be spec- 
ularly reflected. This peculiar property is a result of the existence 
of two bands (conduction and valence) close to the Fermi energy in 
graphene. In the low-energy transport regime, i.e. quasiparticle ener- 
gies E of order 0{A), retroreflection dominates if Ef 3> A, while 
specular reflection dominates if Ef <ti A. 
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re < 1 




re > 1 








FIG. 5: (Color online) The effective impedance caused by FVM il- 
lustrated schematically for all possible cases of smaller, equal, and 
larger Fermi velocity in the normal part of the system. Except for the 
case when the Fermi velocities are identical in the normal and super- 
conducting part of the system (re — 1), part of the Fermi surface does 
not participate in the scattering processes, resulting in a reduction in 
conductance. For instance, when re < 1, total reflection occurs at 
angles ^jv > asin(re), such that FVM effectively acts as a source of 
normal reflection^'*. 



B. Scattering processes 

Consider the case of zero external magnetic field. The 
full Bogoliubov-de Gennes (BdG) equation for the 2D sheet 
graphene normal/s-wave superconductor junction in the xy- 
plane then reads^-^ 



H -Epl 
, All 




= E 



(3) 



where E is the excitation energy, and {u, u} denoting the 
electron-like and hole-like exictations described by the wave- 
function. Assuming that the superconducting region is located 
at a; > and neglecting the decay of the order parameter in the 
vicinity of the interface^^, we may write for the spin-singlet 
order parameter 



Ak = A(0)e"'e(x), 



(4) 



where 8(x) is the Heaviside step function, d is the phase 
corresponding the globally broken U{1) symmetry in the su- 
perconductor, while 6 = atan(fcj,/fc^) is the angle on the 
Fermi surface in reciprocal space (we have adopted the weak- 
coupling approximation with k fixed on the Fermi surface). 
Note that in contrast to previous work, we allow for the pos- 
sibility of unconventional superconductivity in the graphene 



layer since Ak now may be anisotropic. We have applied 
weak-coupling limit, the momentum k is fixed on the Fermi 
surface, such that Ak only has an angular dependence. Since 
we employ a spin-singlet even parity order parameter, the con- 
dition A (6') = A(7r + 6) must be fulfilled. The single-particle 
Hamiltonian is given by 



H 




H± = ~ivp{(7xdx ± (Tydy). (5) 



Here, v-p is the energy-independent Fermi velocity for 
graphene, while (Ji denotes the Pauli matrices. For later use, 
we also define the Pauli matrix vector a — {ax,(Jy,<Jz)- 
These Pauli matrices operate on the sublattice space of the 
honeycomb structure, corresponding to the A and B atoms, 
while the ± sign refers to the two so-called valleys of K and 
K' in the Brillouin zone. The Dirac points earn their sobri- 
quet as valleys from the geometrical resemblance of the band 
dispersion to the aforementioned. The spin indices may be 
suppressed since the Hamiltonian is time-reversal invariant. 
In addition to the spin degeneracy, there is also a valley de- 
generacy, which effectively allows one to consider either the 
one of the H± set. Therefore, the 8 matrix BdG-equation Eq. 
(|3]l reduces to a 4 x 4 matrix BdG-equation, namely 



H± - Epl Akl 

Epl-H± 




E 



(6) 



where have explicitly used that TH = HT . Let us then con- 
sider such that one may write 



p • o- - Epl Akl 



Epl — p • o- 




= E 



(7) 



In the above Hamiltonian, we have only included diagonal 
terms in the gap matrix, i.e. Ak = Akl. This corresponds 
to exclusively intraband-pairing on each of the sublattices A 
and B. In recent work by Black-Schaffer and Doniach^^, it 
was shown that by postulating interband spin-singlet hopping 
between the sublattices, one could achieve dominant d-wave 
pairing in intrinsic graphene. While an onsite attractive po- 
tential is sufficient to achieve s-wave pairing, leading to a di- 
agonal gap matrix, nearest-neighbor interactions couples the 
two sublattices and should yield off-diagonal elements in the 
gap matrix. In this work, we restrict ourselves to anisotropic 
superconducting pairing with diagonal elements in the gap 
matrix, although one would have to take into account off- 
diagonal elements as well for a completely general treatment. 
We comment more on this later. 

Consider an incident electron from the normal side of the 
junction {x < 0) with energy E. For positive excitation en- 
ergies E > 0, the eigenvectors and corresponding momentum 
of the particles read 

^Pl = [l,e'^0,0]Vf'^°^^^ = {E + Ep)/v,, (8) 

for a right-moving electron at angle of incidence 9 (see Fig. |6] 
while a left-moving electron is described by the substitution 
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6 ^ TT — 9. If Andreev-reflection takes place, a left-moving 
hole is generated with an energy E, angle of reflection 6 a, and 
corresponding wave-function 

V-^L = [0,0,l,e-'**]^e-'P'''=°'^^*^ = iE~EF)/vF, (9) 

where the superscript e (h) denotes an electron-Uke (hole- 
like) excitation. Since translational invariance in the in- 
direction holds, the corresponding component of momentum 
is conserved. This condition allows for determination of the 
Andreev-reflection angle 6\ through p*" sin 0a = p'^sin6'. 
From this equation, one infers that there is no Andreev- 
reflection (9a — ±7r/2), and consequently no subgap con- 
ductance, for angles of incidence above the critical angle 

\E- E, 



On the superconducting side of the system (x > w), the pos- 
sible wavefunctions for transmission of a right-moving quasi- 
particle with a given excitation energy E > reads 

(u{9+),u{9+)e'^^,v{e+)e-"^^,v{9+) 



^iq cos p ' X 



X e 



(v{9-), t>(0-)e"^" , u{9-)e-"^' , u(0-)e'(''"""*")' ^ 
, = {E^ - ^E^-\A{9-)\^)/v,. 



iq cos u X 



(11) 



The coherence factors are, as usual, given byai 



u{9) 



v{9) 





^E^ 


-\Ai9W 


E 




-\Ai9W 




E 


= t 


)|, 9- = 



(12) 



9^, and 



ei0 ^ e"^A{9^)/\A{9^)\. The transmission angles 6i^'^ for 
the electron-like (ELQ) and hole-like (HLQ) quasiparticles 
are given by g^'' sin 6*5^ — sin 9, i=e,h. Note that for subgap 
energies E < A, there is a small imaginary contribution to 
the wavevector, which leads to exponentional damping of the 
wavefunctions inside the superconductor. The physical reason 
for this is that there can be no transmission of quasiparticles 
into the superconductor for subgap energies. Note that for 
all wavefunctions listed in the equations above, we have for 
clarity not included a common phase factor e'*^"^ which cor- 
responds to the conserved momentum in the y-direction. A 
possible Fermi vector mismatch (FVM) between the normal 
and superconducting region is accounted for by allowing for 
Ep Ep. The case i?p 3> A corresponds to a heavily doped 
superconducting region. It is also straight-forward to obtain 
the eigenfunctions for the case when the gap matrix consists 
of off-diagonal elements, as opposed to the gap matrix treated 
here with diagonal elements. In particular, for 



(13) 




the eigenfunctions may be obtained from Eq. (fT4b simply by 
switching the phase-factors as follows: 



X e'«'™''«^% ^ {El. + ^E^ -\A{9+)\^)/vp, 

= (^v{9-),v{9-)e'''' ,u{9~)e''^'^'^'^'\u{9-)e'"f''y 



X e 



iq cos 9 X h / rpl 



q"^ = {E'p - ^E^-\A{9-)\^)/ 



Vp. 



(14) 



At this stage, it is appropriate to insert the restriction which 
will be used throughout the rest of this paper, namely A <C 
-Ep. Since we are using a mean-field approach to describe the 
superconducting part of the Hamiltonian, it is implicitly un- 
derstood that phase-fluctuations of the order parameter must 
be smalP''. For this criteria to be fulfilled, the superconducting 
coherence length ^ must be large compared to some character- 
istic length scale of the system^^. Following Ref.'28', the criti- 
cal temperature Tk at which long-range phase-fluctuations of 
the order parameter destroys the ordering when approaching 
the critical temperature from below, is given by 



TK = rf^(l-|r 



(15) 



where |r| ^ , D \s the dimensionality of the system, and 
is the critical temperature predicted by mean-field theory. 
(For an extensive treatment of the effect of phase fluctuations 
in extreme type-II superconductors, see Ref. 29.) Thus, only 
for |r| ^ 1 mean-field theory is a viable option for describ- 
ing superconductivity in the system, corresponding to a large 
coherence length ^. Notice that the Ginzburg temperature Tg, 
which describes the regime where amplitude-fluctuations of 
the order parameter become important (T > To), satisfies 
jtMF > > 2k a natural choice of characteristic length 
scale for the system in the normal state is obviously the Fermi 
wavelength Ap — 2TTVp / Ep, such that the criteria for valid- 
ity of mean-field theory reads ^/Ap ^ 1, or equivalently, 
E^ > A. 

The relevant scattering processes at the N/S graphene in- 
terface are shown in Fig. |6] in the two cases of zero barrier 
and an insulating interface of width w. In the former case, the 
boundary conditions dictate that iIj\x=o — '^\x=o, where 



(16) 



while in the latter case, one must match the wavefunctions at 
both interfaces: 



tp\x=0 = '0lU=O, '4'l\x=w = ^'sU=t, 



(17) 



where we have defined the wavefunction in the insulating re- 
gion 



(18) 



The wavefunctions tp differ from ip in that the Fermi energy 
is greatly shifted by means of e.g. an external potential, such 
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that Ep ^ — Vo where Vq models the potential barrier 
(equivalent to the role of Z in Ref. ,33,) . Also, note that the 
trajectories of the quasiparticles in the insulating region, de- 
fined by the angles 9 and 9a, differ by the same substitution, 
meaning 

sine/ sing = {E + Ef)/{E + E^ - Vq), 
sin 0a/ sin 61 = {E + Ep)/ {E - Ep + Vq). 

(19) 

Finally, note that the subscript ± on the wavefunctions in the 
normal region indicates the direction of their group velocity, 
which in general is different from the direction of momentum, 
as discussed previously. Consequently, although the Andreev- 
reflected hole wavefunction carries a subscript "— " above, one 
should keep in mind that for normal Andreev reflection, the 
direction of momentum is opposite to the group velocity for 
the hole. 
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(T) Incident conduction-band electron (4) Specularly AR mlence-band hole 

(2) Normally AR conduction-band hole (5) Transmitted ELQ 

(3) Reflected conduction-band electron (6) Transmitted HLQ 

FIG. 6: (Color online) The scattering processes taking place at an 
N/S or N/I/S graphene junction. In the former case, the insulating 
region is completely absent, and only the six depicted processes take 
place. Note that only normal Andreev reflection or specular Andreev 
reflection takes place at any given energy E, never both. For an N/I/S 
graphene junction, there are transmitted and reflected electrons and 
holes in the insulating region corresponding to ^i, not shown in the 
above figure. 

Before we go on to presenting results, we make one con- 
ceptual remark. The proximity effect means that an otherwise 
normal system becomes superconducting by virtue of having 
the superconducting wavefunction from a nearby supercon- 
ductor leak into the normal system, thus making it supercon- 
ducting in some region. This is a result of a boundary condi- 
tion imposed on the normal system from the proximate host 
superconductor. The resulting wavefunction in the proximity 
region of the normal system is then a BCS type wavefunction. 
Such a wavefunction unquestionably describes a system with 
a gapped Fermi surface (possibly with nodes on the Fermi- 
surface). It matters not by what microscopic mechanism such 
a state was established, as long as it is there. The effective in- 
teraction giving rise to proximity induced superconductivity 
in graphene close to the surface in contact with an intrinsi- 
cally superconducting host system, is obtained by considering 



the complete superconductor-graphene system and integrating 
out the electrons on the superconducting side. The electrons 
in graphene then experience an effective attractive interaction 
Aeff giving rise to a gap by virtue of hopping into and out of the 
superconducting side. We thus have, by such tunneling pro- 
cesses, Aeff 7^ even if A = 0. Here, A is the electron-electron 
coupling constant giving rise to superconductivity in graphene 
per se. Since graphene intrinsically is a normal system, and 
is well approximated by non-interacting electrons, this cou- 
pling constant vanishes, A = 0. The relationship between the 
gap in the normal region A and (//) is thus A = Xes{ff), 
and this gives a nonzero gap in the vicinity of the proximate 
host superconductor Here, / are fermion annhilation opera- 
tors, and (//) thus represents the pair-amplitude induced in 
the normal graphene region. The proximity-induced gap van- 
ishes rapidly as one goes away from the surface and into the 
bulk of graphene, since Aeff vanishes rapidly as we move away 
from the proximate host superconductor It would in principle 
be incorrect to assert that in the proximity -region of graphene, 
we could have a non-zero anomalous Green's function (//), 
but no gap A, by using a self-consistency relation of the type 
A = A(//) with A = 0^^". Such a self-consistency relation 
does not exist in a normal system which does not supercon- 
duct by itself. However, in a situation where the intrinsic 
A = 0, it could well turn out to be the case that Aeff is small, 
leading to a gap which is very small^^. In our paper, we have a 
thin film graphene system with a bulk superconductor in con- 
tact with the film, deposited on top of the film. If the thick- 
ness of the graphene-film is smaller than the coherence length 
of the bulk superconductor, one obtains a proximity-induced 
superconductiving gap throughout the film. As we shall see, 
our results are quite sensitive to the presence of even a small 
induced gap in graphene. 

III. CONDUCTANCE SPECTRA 

In what follows, we describe how the conductance spectra 
of a N/S and N/I/S graphene junction may be obtained. Ac- 
cording to the BTK formalism^, the normalized conductance 
is given by 

1 r''^ / 

G{eV)^— d9cos9(l^\r{eV,0)\'^ 

where r and are the reflection coefficients for nor- 
mal and Andreev reflection"*", respectively, while Gn is a 
renormalization constant corresponding to the N/N metallic 
conductance^^, 

Gn=/ dgcosg Tn'y2 - (21) 

In this case, we have zero intrinsic barrier such that Z = 
0. We will apply the usual approximation \rA{—eV,6)\ — 
\rA{eV,9)\, which may be shown to hold for a quite general 
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parameter regime. For perfect normal reflection (|rp = 1), 
there is no conductance, while for perfect Andreev reflection 
(|rAp = 1), the conductance is doubled compared to the N/N 
case. In order to obtain these coefficients, we make use of the 
boundary conditions described in the previous section. The 
analytical solution and behaviour of the conductance differs 
in the N/S and N/I/S case, and we proceed with a separate 
treatment of these scenarios. 



A. N/S junction 

Solving the boundary conditions for the wavefunctions at 
the interface leads to the analytical expressions for the reflec- 
tion coefficients: 



2cos9[C+v{e+)v{9-)& 



C^u{9+)u{9-)] 



TA = 



v{9+)v{9~)&('t'--'i'+')Y^ - u{9+)u{9-)Y+ 
2e-"^^ cos9[C+u{9-)v{9+) - C-u(6'-)w(6'+)] 



- 1, 



t;(6l+)w(6l-)e'("^'--'>+)y_ - u{9+)u{e-)Y+ 



(22) 



where we have defined the auxiliary quantities (± 



andY± = (^{c" 



). The interplay between the 



different phases felt by the ELQ and HLQ in the supercon- 
ductor in the case of an anisotropic order parameter enters 
above through the factor e'^''^ K It remains, however, to 
be clarified how this interplay manifests itself in the tunneling 
conductance. Before investigating this in more detail, let us 
briefly consider the isotropic s- wave case first, /.e. A{9) — A, 
such that e'('^"-'^'^) = 1. 



1. Conventional s-wave pairing 



to 



For conventional superconducting pairing, Eq. ( |22] ) reduces 



''A 



(23) 



This case was first studied by Beenakker^. For consistency 
and completeness, we reproduce the results of Ref. 2 [see 
Fig. |7^)]. We point out that Eq. ( |23] ) are valid for any pa- 
rameter range, and not restricted to the heavily doped case 
treated in Ref. 2^. To illustrate the difference, we consider the 
regime E!p = shown in Fig. |7j)). In this case, the standard 
situation of perfect Andreev reflection for subgap energies is 
recovered, with a sharp drop at the gap edge corresponding to 
the onset of quasiparticle transmittance into the superconduc- 
tor. 



2. Anisotropic d-wave pairing 

To treat an unconventional superconducting order parame- 
ter, we must revert to the general expressions in Eq. ( l22b . In 
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FIG. 7: (Color online) Conductance spectra for graphene in a) with 
Ef/A — 10^,andinb)with_Ep — E^. The spectra in a) are identical 
to the result of Ref. 2. In b), the subgap conductance is always close 
to 2Gn, but becomes more constant for increasing Ef since 6^ 
■k/2. We have plotted the ratios {50, 100, 1000} of Ev/A in b), 
from bottom to top. 



order to account for the effect of an anisotropic gap, we choose 
the (ij.2_y2-gap from Tab. HI which in the weak-coupling ap- 
proximation (|k| — fcp) reads A{9) = A cos(26' — 2a). Here, 
a models the relative orientation of the gap in k-space with 
respect to the interface normal as illustrated in Fig. [8] 



Interface 









" rfx2-j,2 gap 




Transmitted ELQ 
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FIG. 8: (Color online) Sketch of 2D N/S graphene junction with an 
anisotropic superconductor. The orientation of the gap in k-space is 
modelled by the angle a. 

We now proceed to investigate how the conductance spec- 
tra of a N/S graphene junction change when going from a s- 
wave to a d-wave order parameter in the superconducting part 
of the system. Consider Fig. |9]for the case of heavily doped 
graphene, where the orientation of the gap is such that the con- 
dition for perfect formation of zero energy states in a metallic 
N/S junction is fulfilled, i.e. A{9) = -A(7r - 9). As shown 
by Tanaka and Kashiwaya^", this gives rise to a quasiparticle 
interference between the ELQ and HLQ since they feel dif- 
ferent phases of the pairing potentials due to their different 
trajectories of transmittance into the superconductor This re- 
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suits in a bound surface states with zero energy^^ close to the 
N/S interface. For the N/S graphene junction studied here, the 
explicit barrier potential is zero, while FVM effectively acts 
as as source of normal reflection. From Fig. |9] one may in- 
fer that a peak at zero bias is present in the presence of FVM, 
although the ZBCP does not increase in magnitude with in- 
creasing FVM. We will later study how the presence of an 
intrinsic barrier in the form of a thin, insulating region sepa- 
rating the normal and superconducting part affects the ZBCP. 
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FIG. 10: (Color online) Conductance spectra for a) doped graphene 
in with Ep — 10"* A for a d-wave order parameter and b) undoped 
graphene with Ep — Ep. In both cases, we have set Ep/A — 100 
and investigate how the conductance spectra evolves upon rotating a 
from to 7r/4 in steps of n/20. In a), it is seen that the peak of the 
conductance shifts from eV — A to progressively lower values as 
a increases. This is in contrast to metallic N/S junctions. In b), we 
plot the conductance for a — 7r/4. We find virtually no difference 
between various orientations of the gap in this case, and the spectra 
are quite similar to the metallic N/S case with zero barrier. 



B. N/I/S junction 



FIG. 9: (Color online) Conductance spectra for doped graphene in 
with Ep = 10*A for a d-wave order parameter with orientation an- 
gle Q = 7r/4. A ZBCP is present for large FVM, and becomes 
unobservable narrow for Ep/ A < 10. For a — 0, the conductance 
spectra are essentially identical to those in Fig. [7] 



Next, we plot the conductance spectra for doped graphene 
to see how they evolve upon a rotation of the gap. The be- 
haviour is quite distinct from that encountered in a N/S metal- 
lic junction. From Fig. [TOh). we see that the peak of the con- 
ductance shifts from eV = A to progressively lower values 
as a increases from to 7r/4. In this respect, the conduc- 
tance spectra actually mimicks a lower value of the gap than 
what is the case, if one were to infer the gap magnitude from 
the position of the singularity in the spectra. As of such, for 
a given FVM, determining the magnitude of the gap by the 
usual method of locating the characteristic feature in the con- 
ductance spectra is not as straight-forwards in N/S graphene 
junctions as in the metallic case. Indeed, multiple measure- 
ments with several different interface orientations would in 
general be required to obtain the correct value of the gap. This 
should be a unambigously observable feature in experiments, 
and provides a direct way of testing our theory. Finally, we 
consider graphene in Fig. [TOb ) with Ep = E^ for a d-wave 
order parameter Upon varying a from to 7r/4, there is now 
little distinction between different angles of orientation. The 
conductance spectra are in this case very resemblant to the 
metallic N/S case for zero barrieriS,. 



We now consider the conductance of an N/I/S graphene 
junction, where I denotes an "insulating" (see introduction) 
region modelled by a very large energy potential for the quasi- 
particles. Solving the boundary conditions introduced in Sec. 
im we obtain analytical expression for r and ta, which is all 
that is required in order to calculate the conductance. How- 
ever, these expressions are very large and the reader may 
consult Appendix A for their explicit form. In the follow- 
ing, we will not work exclusively in the thin barrier-limit 
d 0, Vq — > 00 as in Ref. 13. Some aspects of including an 
insulation region of arbitrary width and strength were very re- 
cently discussed in Ref. J4, albeit only in the case of isotropic 
s-wave pairing. We now treat the two cases of s-wave and 
d-wave pairing separately. 



1. Conventional s-wave pairing 

This case was first studied by Bhattacharjee and 
Sengupta'^*. To quantify the parameters in the insulating re- 
gion, we will measure the width w of region I in units of Ap 
and the potential barrier Vq in units of Ep. First, we briefly 
show that we are able to reproduce the qualitative findings of 
Ref. 13. As shown in Appendix A, it is convenient to intro- 
duce the parameter x = Vqw/vp in the thin-barrier limit. In 
this case, the reflection coefficients r and ta exhibit an inter- 
esting oscillating behaviour as a function of x- To see this, 
consider Fig. [TT] where we have plotted the voltage depen- 
dence of the normalized conductance for several values of x- 
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For X = 0, we reproduce the result of Fig. [TJ)). This is reason- 
able since the conductance of an N/l/S junction with x — 0, 
i.e. zero width, should be the same as an N/S junction. The tt- 
periodicity is reflected in that the curves for x = and x — 
are identical. For x rt = 0, 1, 2, .., there is a source of 

normal reflection at the interfaces due to the insulating region, 
and consequently the subgap conductance is reduced from its 
ballistic value 2Gn. Even in the presence of a FVM, Ep ^ Ep, 
the spectra of Fig. [TT| retain their 7r-periodicity. However, the 
FVM acts as a source of normal reflection such that one does 
not have nearly perfect Andreev reflection at subgap energies. 

Our results differ slightly from those reported in Ref. Ii3[ 
Although we obtain qualitatively exactly the same depen- 
dence on X of the conductance, it is seen by comparing our 
Fig. [TT]with Fig. 1 of Ref. 1 3 that our curves are phase shifted 
by 7r/2 in X in comparison. As a consequence, we regain the 
N/S conductance result when x = instead of x = 7r/2 as 
reported in Ref. [H!. Physically, this seems to be more reason- 
able since x = corresponds to the case of an absent barrier, 
a situation where there is no source of normal reflection be- 
sides the condition that momentum in the direction parallell 
to the barrier must be conserved. We have also verified that 
our X = result coincides with the results obtained using 
the full expressions (see Appendix) without assuming a thin- 
barrier limit when we let both w and Vb go to zero. We believe 
that this minor discrepancy between our results and the results 
of Ref. 13 stems from a sign error in their Eq. (5) and also in 
their expression for fcf,(A;^) in the text above Eq. (5). 



reflection that takes place at zero bias voltage, thus reducing 
the conductance. 
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FIG. 12: (Color online) Plot of the zero-bias conductance as a func- 
tion of X with Ep/A = 100 for a N/I/S graphene junctions with 
s-wave pairing. 
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FIG. 1 1 : (Color online) Conductance spectra for an N/I/S graphene 
junction with -Ep/A = Ef/A — 100, using s-wave pairing. We 
reproduce the same results as Ref. 13 with a 7r-periodicity in the pa- 
rameter X- However, we obtain a phase shift of 7r/2 in x compared 
to their results. We believe that this difference pertains to a minor 
sign error in the wavefunctions used in Ref. 

To unveil the periodicity even more clear, consider Fig. [12] 
for a plot of the zero-bias conductance as a function x- The 
cases Ep ~ Ep and Ep ^ Ep display a striking difference. 
The qualitative shape of the curves is equal, but the amplitude 
is diminished with increasing FVM. This may in similarity 
to the above discussion be attributed to the increased normal 



2. Anisotropic d-wave pairing 



We now contrast the s-wave case with an anisotropic pair- 
ing potential to see how the spectra are altered. Consider first 
Fig. [13] for a plot of the tunneling conductance in the un- 
doped case. We consider the two angles a = and a — tt/A 
as representatives for the two types of qualitative behaviour 
that may be expected in a d-wave superconductor/normal 
graphene junction. The latter corresponds to perfect formation 
of ZES in the metallic counterpart junction. From the spectra, 
one infers that for a = 0, tunneling into the nodes of the su- 
perconducting gap destroys the nearly perfect Andreev reflec- 
tion for subgap energies obtained in Fig. [TT] When a = 7r/4, 
one observes the formation of a ZBCP which peaks at twice 
the normal state conductance. It is also interesting to note that 
the zero bias conductance remains unchanged upon increasing 
X- Therefore, the equivalent of Fig. [T2]in the present d-wave 
case is G{0) ~ 2G'n, regardless of x- 

Introducing a FVM between the superconducting and nor- 
mal parts of the system, the spectra are rendered less sensitive 
to any increase in x> as seen in Fig. [14] For a = 0, the 
spectra are essentially identical to the doped s-wave case. For 
a = 7r/4, it is seen that the formation of a ZBCP becomes 
even more protruding, and that the zero bias conductance is 
still insensitive to any increase in x- Therefore, one is led to 
conclude that the normalized zero bias conductance G(0) / Gn 
in the d-wave case is constantly equal to nearly 2, regardless 
of X and the magnitude of the FVM. 
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FIG. 13: (Color online) Conductance spectra for an N/I/S graphene 
junction with E^/A = Ev/ A — 100, using d-wave pairing. 
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FIG. 14: (Color online) Conductance spectra for an N/I/S graphene 
junction with Ep/Ep = 10 and Ep/A = 100, using d-wave pairing. 



IV. CONDUCTANCE-OSCILLATIONS 

In this section, we investigate a feature of the conductance 
spectra that is in common for both the s-wave and d-wave 
case: an oscillatory behaviour as a function of applied bias 
voltage. Consider first a N/I/S graphene junction. In the thin- 
barrier limit defined as w ^ and Vb — > oo with s-wave 
pairing, Ref. J3i reported a 7r-periodicity of the conductance 
with respect to the parameter x = Vqw/vp, as discussed in 
the previous section. We now show that by not restricting 
ourselves to the thin-barrier limit, new physics emerges from 
the presence of a finite-width baiTier. We measure the width w 
of region I in units of Ap and the potential barrier Vq in units of 
Ep. The linear dispersion approximation is valid^' up to ~ 1 
eV, and we will consider Fermi energies in graphene'* ranging 
from the undoped case iJp ~ meV to Ep w 100 meV in the 



doped case, setting the gap value to A = 1 meV. Owing to 
the restriction of Ep :$> A, we fix Ep — 100 A, and also set 
Vq ~ 500A in order to model the effective potential barrier. 

Consider Fig. [15] where we plot the normalized tunneling 
conductance in case of s-wave pairing, for both a doped and 
undoped normal part of the system. The most striking new 
feature compared to the thin-barrier limit is the strong oscilla- 
tions in the conductance as a function of eV. For subgap en- 
ergies, we regain the N/S conductance for undoped graphene 
when X 0, with nearly perfect Andreev reflection. The 
same oscillations are seen in the d-wave pairing case, shown 
in Fig. [16] To model the d-wave pairing, we have used the 
d^2_y2 model A(e') = Acos(26' - 2a) with a = tt/4. The 
parameter a effectively models different orientations of the 
gap in k-space with regard to the interface, and a = 7r/4 
coiTesponds to perfect formation of ZES in N/S metallic junc- 
tions. For a = 0, the d-wave spectra are essentially identical 
to the s-wave case, since the condition for formation of ZES 
is not fulfilled in this case^". It is seen that in all cases shown 
in Figs. [15] and [16] the conductance exhibits a novel oscilla- 
tory behavior as a function of applied bias voltage eV as the 
width w of the insulating region becomes much larger than the 
Fermi wavelength, i.e. w ^ Xp. 

The oscillatory behavior of the conductance may be under- 
stood as follows. Non-relativistic free electrons with energy 
E impinging upon a potential baiTier Vq are described by an 
expontentially decreasing non-oscillatory wavefunction e'*^^ 
inside the baiTier region if E < Vq, since the dispersion es- 
sentially is fc ~ y/E — Vq. Relativistic free electrons, on the 
other hand, have a dispersion k ~ (i? — Vo), such that the 
corresponding wavefunctions do not decay inside the barrier 
region. Instead, the transmittance of the junction will display 
an oscillatory behavior as a function of the energy of inci- 
dence E. In general, a kinetic energy given by ~ /c" will 
lead to a complex momentum k ^ (E — Vb)^/" inside the 
tunneling region, and hence damped oscillatory behavior of 
the wave function. Relativistic massless fermions are unique 
in the sense that only in this case (a = 1) is the momentum 
purely real. Hence, the undamped oscillatory behavior at sub- 
gap energies appears as a direct manifestation of the relativis- 
tic low-energy Dirac fermions in the problem. This observa- 
tion is also linked to the so-called Klein paradox which occurs 
for electrons with such a relativistic dispersion relation, which 
has been theoretically studied in normal graphene^^. 

We next discuss why the illustrated conductance spectra 
are different for s-wave and d-wave symmetry, in addition to 
comparing the doped and undoped case. The difference in 
doping level between the superconducting and normal part of 
the system may be considered as an effective FVM, acting as 
a source of normal reflection in the scattering processes. This 
is why the subgap conductance at thin barrier limit is reduced 
when Ep ^ Ep. Moving away from the thin barrier limit, it is 
seen that oscillations emerge in the conductance spectra. For 
s-wave pairing, the amplitude of the oscillations is larger for 
Ep ^ Ep than for the case of no FVM, but the period of os- 
cillations remains the same. This period depends on w, while 
the amplitude of the oscillations is governed by the wavevec- 
tors in the regions I and S. The maximum value of the oscil- 



12 



lations occurs when 2w equals an integer number of wave- 
lengths, corresponding to a constructive interference between 
the scattered waves. Physically, the amplitude-dependence of 
the oscillations on the doping level originates with the fact that 
any FVM effectively acts as an increase in barrier strength. By 
making Vq larger, one introduces a stronger source of normal 
reflection. When the resonance condition for the oscillations 
is not met, the barrier reflects the incoming particles more ef- 
ficiently. This is also the reason why increasing Vq directly 
and increasing the FVM has the same effect on the spectra. 

We now turn to the difference between the s-wave and d- 
wave symmetries. It is seen that the conductance is reduced 
in the d-wave case compared to the s-wave case, and is ac- 
tually nearly constant for Ep — 0. One may understand 
the reduction in subgap conductance in the undoped case as 
a consequence of tunneling into the nodes of the gap, which 
is not present in the s-wave case. Hence, Andreev reflection 
which significantly contributes to the conductance, is reduced 
in the d-wave case compared to the s-wave case. Moreover, 
we see that a ZBCP is formed when Ep ^ Ep, equivalent to a 
stronger barrier, and this is interpreted as the usual formation 
of ZES leading to a transmission at zero bias with a sharp drop 
for increasing voltage. 
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FIG. 16: (Color online) Same as Fig. [15] but now for d-wave pairing. 
We have fixed Vo /A = 500 and E'p/A = 100. 



V. DISCUSSION 
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FIG. 15: (Color online) Tunneling conductance of N/I/S graphene 
junction for s-wave pairing in the undoped and doped case (see main 
text for parameter values). We have fixed Vo/ A = 500 and iJp/ A = 
100. It is seen that for increasing w, a novel oscillatory behaviour of 
the conductance as a function of voltage is present in all cases. 



By means of the proximity effect, Heersche et al. suc- 
cessfully induced superconductivity in a graphene layer^ (see 
also Ref. [Sj). This achievement opens up a vista plethora 
of new, exciting physics due to the combination of the pe- 
culiar electronic features of graphene and the many inter- 
esting properties of superconductivity. For our theory to be 
properly tested experimentally, it is necessary to create N/S 
and N/I/S graphene junctions. Junctions involving normal 
graphene with insulating regions have recently been experi- 
mentally realized'*-^. In our work, we have discussed novel 
conductance-oscillations in a N/l/S graphene junction that 
arise when moving away from the thin-barrier limit discussed 
in Ref. 13. While reaching the thin-barrier limit might pose 
some difficulties from an experimental point of view, our 
predictions are manifested when using wide barriers, which 
should be technically easier to realize. In order to reach the 
doped regime, this could be achieved by either chemical dop- 
ing or using a gate voltage to raise the Fermi level in the super- 
conducting region^^i^. The relevant magnitudes for the vari- 
ous physical quantities present in such an experimental setup 
has been discussed in the main text of this paper. 

It is also worth mentioning that since we have assumed a 
homogeneous chemical potential in each of the normal, in- 
sulating, and superconducting regions, the experimental re- 
alization of the predicted effects require charge homogeneity 
of the graphene samples. This is a challenging criteria, since 
electron-hole puddles in graphene imaged by scanning single 
electron transistor^^ suggest that such charge inhomogeneities 
probably play an important role in limiting the transport char- 
acteristics of graphene-"^. In addition, we have neglected the 
spatial variation of the superconducting gap near the N/S in- 
terfaces. The suppression of the order parameter is expected 



to least pronounced when there is a large FVM between the 
two regions. However, the qualitative results presented in this 
work are most hkely unaffected by taken into account the re- 
duction of the gap near the interface. 

VI. SUMMARY 

In summary, we have studied coherent quantum 
transport in normal/superconductor (N/S) and nor- 
mal/insulator/superconductor (N/I/S) graphene junctions, 
investigating also the role of d-wave pairing symmetry on 
the tunneling conductance. We elaborate on the results 
obtained in Ref. 18, namely a new oscillatory behaviour of 
the conductance as a function of bias voltage for insulating 
regions that satisfy w ^ Ap, which is present both for s- and 
d-wave pairing. This is a unique manifestion of the Dirac-like 
fermions in the problem. In the d-wave case, we have studied 
the conductance of an N/S and N/I/S junction in order to 
make predictions of what could be expected in experiments, 
providing both analytical and numerical results. We find very 
distinct behaviour from metallic N/S junctions in the presence 
of a FVM: a rotation of a is accompanied by a progressive 
shift of the peak in the conductance, without any formation of 
a ZBCP except for a — tt/4. All of our predictions should be 
easily experimentally observable, which constitutes a direct 
way of testing our theory. 
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APPENDIX A: NORMAL- AND ANDREEV-REFLECTION 
COEFFICIENTS FOR N/VS JUNCTIONS 

Solving the boundary conditions Eq. ( [TtT i. we obtain the 
following expressions for the normal reflection coefficient r 
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and the Andreev-reflection coefficient r^: 

r = U{A + C) + th{B + D)-l, 
rA = UA' + C') + th{B' + D'), (Al) 
where the transmission coefficients read 

te = 2cos6i[c-"'-*(B' + D') - {B'c-'^^ - D'e'^'')]p-\ 
th = te[e"''*(^'e-'«^^ - C'e'^^) - A' - C] 

X [B' + D' ^ e'^-^(S'e^'^~^ - D'e'^^)]^\ (A2) 

with the definition 

p = [e-'^-'iB' + D') - (B'e-'^^ - D'e'^-")] 
X [e-'^iA + C) + (Ac'^~ - Ce-^^)] 
- [{De-'^ - Be'^) - e''%B + D)] 
X [{A'e-'~^^ - C'e'^^ ) - e""'^ {A' + C')] (A3) 

We have defined the auxiliary quantities 

A = u+e'('?^~f^)[l- (e'^^-e'^^)(2cos^)"i], 
B = t;_e'(«"-P^)[l- (e'^~-e'^")(2cos^)^i], 
C = u+e'(f^+«^)(e'^" - e'^^ ) (2 cos 6i)-\ 
D = t;_e'(P^+«")(e'^^ - e'^")(2cos6i)-\ (A4) 

and similarly introduced 

A' = w+e'^^^+f-^^^il + (e'^"" - e-'*~^)(2cos^^)-i], 
B' = u_e'('?"+f"-'^")[l + (e'^" - e-'^^-^)(2cos^A)"^], 
C = v+e'^'i^-P'"l'*\e-'^^ cose 
D' = u-c'^'^'-P'-''''\e-'^^ - e"'")(2cos^A)"^ (A5) 

For more compact notation, we have finally defined 

= q'^ cos 9^w, q^ ~ q^ cos 6~w, 

= p'^ cos9w, p~ = p^ COS0AW. (A6) 

In the thin-barrier limit defined as w ^ and Vq oo, one 
may set 

9-^0, 9a^0, q±^ 0, p± ^ TX, (A7) 
where x = Vow/vf. 
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not affect the results of Ref. [T^ in any manner since the case 
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